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Abstract: 

At high temperatures, strongly interacting matter becomes a plasma of deconhned quarks 
and gluons. In statistical QCD, deconfinement and the properties of the resulting quark- 
gluon plasma can be investigated by studying the in-medium behaviour of heavy quark 
bound states. In high energy nuclear interactions, quarkonia probe different aspects of the 
medium formed in the collision. We survey the results of recent charmonium production 
studies in SPS and RHIC experiments. 
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1. Introduction 


The states of matter, their defining features and the transitions between them have always 
been among the most challenging problems of physics. Strongly interacting matter opens 
a new chapter for such studies. The elementary particles of strong interaction physics, the 
hadrons, make up most of the observable matter in the universe. They have intrinsic sizes 
and masses, and the ultimate aim of the theory of strong interactions is to derive these. 
Quantum chromodynamics (QCD) does that, in terms of massless gluons and almost 
massless quarks - it is their dynamics which generates the scales. Thus the nucleon is 
a bound state of three very light quarks; the kinetic and potential energy of the system 
provide the nucleon mass, and the binding potential determines the nucleon radius. 

The interaction between quarks is based on their intrinsic colour charge, just as that 
between electrons and protons or nuclei is determined by their electric charge. The form 
of the interaction is quite different, however. The Coulomb potential vanishes for large 
separation distance, so that electric charges can be separated and have an independent 
existence. In contrast, the potential between quarks increases with separation, so that an 
infinite energy would be needed to isolate a quark. In other words, the quark constituents 
of a hadron are confined and not just bound. 

Strong interaction thermodynamics shows that quark confinement has its inherent limits. 
As bound states of quarks, hadrons have an intrinsic size, the binding radius. At extreme 
density, when several hadrons are compressed into a spatial volume normally occupied 
by a single hadron, it becomes impossible to identify among all the overlapping states 
a specific quark-antiquark pair or quark triplet as belonging to a certain hadron. The 
medium becomes instead a dense multi-quark environment in which any quark can move 
as far as it wants, since it is always close enough to other quarks to satisfy confinement 
conditions. Confinement as a large distance feature thus loses its meaning in the short 
distance world of matter at extreme density. 

It is therefore expected that with increasing temperature, strongly interacting matter 
will undergo a transition from a hadronic phase, in which the constituents are colour- 
neutral bound states, to a plasma of deconfined colour-charged quarks and gluons. This 
deconfinement transition in QCD is quite similar to the insulator-conductor transition in 
normal condensed matter physics: it corresponds to the onset of colour conductivity. 

The thermodynamics of strongly interacting matter can be derived through first principle 
calculations in finite temperature lattice QCD, and we shall return to these studies in the 
next section. They provide quite detailed information on the temperature dependence of 
thermodynamic observables, such as pressure or energy density, and they show how these 
quantities change as the system passes through the deconfinement transition. This raises 
a crucial question: if we were given a box containing strongly interacting matter, how 
could we determine the temperature of the medium and how could we specify the state 
in which this medium is in? This question is the central topic of our survey, and we will 
show that the in-medium behaviour of heavy quark bound states provides an excellent 
tool to answer it. 

Since some thirty years it is known that besides the almost massless up and down quarks 
of the everyday world, and the still relatively light strange quarks required to account for 
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the strange mesons and hyperons observed in hadron collisions, there are heavy quarks at 
the other end of the scale, whose bare masses alone are larger than those of most of the 
normal hadrons. These heavy quarks first showed up in the discovery of the J/ip meson 
|T], of mass of 3.1 GeV; it is a bound state of a charm quark (c) and its antiquark (c), 
each having a mass of some 1.2-1.5 GeV. On the next level there is the T meson [2J, with 
a mass of about 9.5 GeV, made up of a bottom or beauty quark-antiquark pair (bb), with 
each quark here having a mass around 4.5 - 4.8 GeV. Both charm and bottom quarks 
can of course also bind with normal light quarks, giving rise to open charm (D) and open 
beauty ( B) mesons. The lightest of these ‘light-heavy’ mesons have masses of about 1.9 
GeV and 5.3 GeV, respectively. 

The bound states of a heavy quark Q and its antiquark Q are generally referred to as 
quarkonia. Besides the initially discovered vector ground states J/ip and T, both the cc 
and the bb systems give rise to a number of other stable bound states of different quantum 
numbers. They are stable in the sense that their mass is less than that of two light-heavy 
mesons, so that strong decays are forbidden. The measured stable charmonium spectrum 
contains the IS' scalar rj c and vector J/ip, three IP states \c (scalar, vector and tensor), 
and the 2 S vector state ip', whose mass is just below the open charm threshold. There 
are further charmonium states above the ip 1 ] these can decay into DD pairs, and we shall 
here restrict our considerations only to quarkonia stable under strong interactions. The 
observed stable charmonium states are summarized in table 1 and the corresponding bot- 
tomonium (b b) states in table 2. The binding energies A E listed there are the differences 
between the quarkonium masses and the open charm or beauty threshold, respectively. 


state 

Vc 

JH> 

XcO 

Xcl 

Xc2 

ip' 

mass [GeV] 

2.98 

3.10 

3.42 

3.51 

3.56 

3.69 

A E [GeV] 

0.75 

0.64 

0.32 

0.22 

0.18 

0.05 


Table 1: Charmonium states and binding energies 


state 

T 

XbO 

Xbl 

Xb2 

'Y' / 

Xbo 

Xbi 

Xb2 

Jff 

mass [GeV] 

9.46 

9.86 

9.89 

9.91 

10.02 

10.23 

10.26 

10.27 

10.36 

A E [GeV] 

1.10 

0.70 

0.67 

0.64 

0.53 

0.34 

0.30 

0.29 

0.20 


Table 2: Bottomonium states and binding energies 


Quarkonia are rather unusual hadrons. The masses of the light hadrons, in particular 
those of the non-strange mesons and baryons, arise almost entirely from the interaction 
energy of their nearly massless quark constituents. In contrast, the quarkonium masses are 
largely determined by the bare charm and bottom quark masses. These large quark masses 
allow a very straightforward calculation of many basic quarkonium properties, using non- 
relativistic potential theory. It is found that, in particular, the ground states and the lower 
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excitation levels of quarkonia are very much smaller than the normal hadrons, and that 
they are very tightly bound. Now deconfinement is a matter of scales: when the separation 
between normal hadrons becomes much less than the size of these hadrons, they melt to 
form the quark-gluon plasma. What happens at this point to the much smaller quarkonia? 
When do they become dissociated? That is the central theme of this report seen from 
another angle. We shall show that the disappearance of specific quarkonia signals the 
presence of a deconfined medium of a specific temperature |3i. Thus the study of the 
quarkonium spectrum in a given medium is akin to the spectral analysis of stellar media, 
where the presence or absence of specific excitation lines allows a determination of the 
temperature of the stellar interior. 

The prediction of a new state of strongly interacting matter led quite naturally to the 
question of where and how to observe it. The very early universe is an obvious case: 
according to the usual evolution equations, giving its energy density as function of age, it 
must have been a quark-gluon plasma in the first 10 ps after the big bang. The core of 
neutron stars is another possibility, and they may in fact consist of deconfined quarks. But 
the essential boost for the field came from the idea to use high energy nuclear collisions to 
form small and short-lived bubbles of the quark-gluon plasma in the laboratory, making 
them accessible to terrestial experimental study. This idea clearly contains assumptions 
to be verified, and we shall return to its basis later on. Here we only note that again 
an absolutely crucial question is how to probe the thermal conditions of the medium 
formed in such collisions. In view of the evolution of such a medium, this is a much 
more complex question than that addressed in the case of a box containing equilibrated 
strongly interacting matter. Nevertheless, as we shall see, quarkonia may well provide a 
tool for this more general task as well. 

After this introduction, we shall first review the essential aspects of deconfinement (section 
2.) and then summarize the basic properties of quarkonia and the dynamics of their disso¬ 
ciation (section 3.). Following this, we survey different approaches to quarkonium binding 
in QCD thermodynamics (section 4.). We then review the main theoretical features of 
quarkonium production in elementary as well as p —A and nucleus-nucleus collisions (sec¬ 
tion 5.), and finally we try to see what one can learn from the present experimental results 
(section 6.). 

2. Colour Deconfinement 


The transition from hadronic matter to a plasma of deconfined quarks and gluons has 
been studied extensively in finite temperature lattice QCD. We shall consider here the 
case of vanishing baryon-number density (baryons and antibaryons in equal numbers) and 
comment only briefly on the situation at finite baryon density at the end of this section. 
At the transition point, the energy density of the system increases by the latent heat 
of deconfinement, i.e., it grows from a value determined by the degrees of freedom of a 
hadron gas to a much higher one governed by the degrees of freedom of a quark-gluon 
plasma. To illustrate this, we recall that the energy density of an ideal gas of massless 
pions is 



( 1 ) 
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corresponding to the three possible pion charges, while an ideal qnark-gluon plasma, with 
two massless quark flavours, gives 

7 2 2 

e QGP = {2x8+ - [2x2x2x3]}^T 4 = 37^T 4 ~ 12 T 4 , (2) 

as determined by the 16 gluonic and 24 quark-antiquark degrees of freedom. Hence in the 
vicinity of T = T c , the energy density increases by more than a factor ten. The actual 
behaviour, as obtained in finite temperature lattice QCD for different flavour compositions 
0EI, is shown in Fig. |T| 




Figure 1: Energy density vs. T [T 5j Figure 2: A = (e — 3 P)/T 4 vs. T flj 

Above T c , the medium consists of deconfined quarks and gluons. We emphasize that 
deconfinement does not imply the absence of interaction - it is only the requirement to 
form colour neutral bound states that has been removed. In Fig.EJit is seen that up to T ~ 
3 T c the ‘interaction measure’ A = (e — 3 P)/T l remains sizeable and does not vanish p], 
as it would for an ideal gas of massless constituents. This strong interaction in the quark- 
gluon plasma has been interpreted in various ways, as gradual onset of deconfinement 
starting from high momenta jjj, as difference between physical and perturbative vacua |S], 
and even in terms of a large number of coloured “resonance” states PJ. Recent studies of 
the relevant degrees of freedom in the quark-gluon plasma nm seem to indicate, however, 
that a description based on quarks having some kind of thermal mass is more likely to be 
correct. 

To specify the nature of the critical behaviour observed at the transition, recall that in 
the limit of large quark mass (rn q —> oo) we recover pure gauge theory and the average 
value of the Polyakov loop, 

L(T) ~ exp{—F^/T}, (3) 

serves as order parameter EU Here 

F QQ = r !™ f qq( t ’ T ) ( 4 ) 

denotes the large distance limit of the free energy of a heavy quark-antiquark pair, which 
diverges in the confined regime and becomes finite through colour screening in the quark- 
gluon plasma, 

f = 0 T <Tl confinement 

L(T) (5) 

[ 7^ 0 T > Tl deconfinement 
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This defines a critical deconfinement temperature T L , separating the confined and decon- 
fined phases. The underlying mechanism for such behaviour is a global Z 3 invariance of 
the Lagrangian of pure SU( 3) gauge theory, which becomes spontaneously broken in the 
deconfined phase, for T > T L . 

In the other extreme, the chiral limit of massless quarks ( m q —> 0), the QCD Lagrangian 
is chirally symmetric, and the chiral condensate 

X(T) = ($$) ( 6 ) 


serves as order parameter. At low temperatures, the massless quarks acquire an effective 
mass by dressing with gluons, so that the chiral symmetry is spontaneously broken. At 
high temperatures, thermal motion removes the dressing, so that the behaviour of y(T), 


X(T) 



T <T X chiral symmetry broken 
T >T X chiral symmetry restored, 


(7) 


defines a chiral restoration temperature T x separating the phases of broken and restored 
chiral symmetry. 

While there thus exists well-defined thermal critical behaviour in the two limits of large 
and small quark mass, the QCD Lagrangian has for general finite quark masses neither 
a global Z 3 (deconfinement) nor a chiral symmetry. Nevertheless it is found that both 
the Polyakov loop and the chiral condensate continue to vary rapidly in the same narrow 
temperature interval, as shown in Fig. |21 for the case of two relatively light quark flavours 
|12j . We can thus still identify an “almost critical” behaviour, the so-called “rapid cross¬ 
over”. Since this occurs for L(T ) and y(T) at the same temperature, deconfinement and 
chiral symmetry restoration are said to coincide. 
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Figure 3: Polyakov loop and chiral condensate vs. temperature [T2| 

The “critical temperature” determined in this way depends on the masses and the number 
of quark flavours. For two light as well as for two light plus one heavier quark flavour, most 
studies D3 indicate T c = 175 ± 10 MeV, as shown in Fig. EJ although some indications 
for a higher value were recently reported m 
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Figure 4: Critical temperature results from different lattice studies [12 j 


The precise form of critical behaviour in the limiting cases also depends on the number 
of colours and flavours. The situation isillustrated in Fig. [21 for colour SU(3) and three 
quark flavours (w, d, s ). In the gauge theory limit rn q —> oo for all flavours, the transition 
is discontinuous (first order). Decreasing the quark mass decreases the discontinuity, 
until it vanishes on a line of second order transitions, beyond which there is the rapid 
cross-over region. For three massless flavours, the chiral transition is also first order, 
and increasing the quark mass leads again to a decreasing discontinuity and a line of 
second order transitions. For two massless quark flavours and a sufficiently massive third 
flavour, the transition is of second order up to the two-flavour limit. The “physical 
point” corresponding to two light u and d quarks and one heavier s quark has now been 
established with some certainty and appears to fall into the cross-over region. 


CO CD 



Figure 5: Critical behaviour as function of quark masses 
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The phase structure at non-vanishing baryon density has stimulated much interest m, 
but it encounters calculational problems in the computer simulation. Quite recently, a 
number of methods to overcome these have been put forward 03. and today we be¬ 
lieve that the phase diagram in T, /j has the form shown in Fig. El with /j, denoting the 
baryochemical potential. 



Figure 6: Critical behaviour for two light quark flavours, as function of T and /j 

We can thus conclude this section by noting that the thermodynamics of strongly inter¬ 
acting matter is today known in considerable quantitative detail, with deconfinement for 
H = 0 setting in at T c = 175 ± 10 MeV. The actual nature of the “transition” remains, 
however, quite enigmatic, since it occurs as “rapid cross-over” and not as a genuine ther¬ 
mal phase transition. Further clarification could come from studies in the context of 
geometric critical behaviour, relating the onset of deconfinement to cluster formation and 
percolation. 

Hence we know that confined hadronic and deconhned quark-gluon states of matter exist. 
How can we probe strongly interacting matter through specific observables, in order to 
determine its thermal parameters (temperature, energy density) as well as to identify 
its confinement status? Since we want to show that this can be done with the help of 
quarkonium states, we first summarize in the next section the essentials of quarkonium 
binding and break-up. 

3. Quarkonium Binding and Dissociation 

We had defined quarkonia as bound states of heavy quarks which are stable under strong 
decay, i.e., M c5 < 2 Mp for charmonia and M bb < 2 Mb for bottomonia. Since the quarks 
are heavy, with m c ~ 1.2 — 1.5 GeV for the charm and m b — 4.5 — 4.8 GeV for the bottom 
quark, quarkonium spectroscopy can be studied quite well in non-relativistic potential 
theory. The simplest (“Cornell”) confining potential mi for a QQ at separation distance 
r has the form 

V(r) = a r — — (8) 

r 

with a string tension a ~ 0.2 GeV 2 and a Coulomb-like term with a gauge coupling 




a ~ 7r/12. The corresponding Schrodinger equation 


( 2mc - —V 2 + V(r)l $,(r) = M^(r) (9) 

l m c ) 

then determines the bound state masses Mj and the wave functions $j(r), and with 

( r i) = J d 3 r r 2 |<hj(r)| 2 / J d 3 r | < f> i (r)| 2 . (10) 

the latter in turn provide the (squared) average bound state radii. 

To obtain a first idea of the results, we consider the semi-classical solution obtained by 
making use of the uncertainty relation ( p 2 )(r 2 ) ~ 1 in the expression for the QQ energy 

n 2 1 

E(r ) = 2m + — + V(r) ~ 2m H-- + V[r). (11) 

m mr 1 

Minimizing E{r ) with respect to r, dE/dr = 0, gives 

r 0 — 0.44 fm M 0 = E(r 0 ) ~ 3.1 GeV (12) 

for the generic charmonium ground state, and 

r 0 ~ 0.33 fm M 0 = E(r 0 ) ~ 9.6 GeV (13) 

and for the corresponding bottomonium state. Here ro specifies the QQ separation dis¬ 
tances, i.e., it gives twice the radius. From Tables 1 and 2 we see that these first estimates 

already reproduce quite well the spin-averaged ground states. 

The exact solution of eq. © gives in fact a very good account of the full (spin-averaged) 
quarkonium spectroscopy, as seen in Table 3p|. The line labelled AM shows the differ¬ 
ences between the experimental and the calculated values; they are in all cases less than 
1 %. Again ro gives the QQ separation for the state in question. The input parameters 
for these results are m c = 1.25 GeV, m b = 4.65 GeV, \fcr = 0.445 GeV, a = 7t/12. 


state 

J/4 

Xc 

if. y 

T 

Xb 

'Y'/ 

x' b 


mass [GeV] 

3.10 

3.53 

3.68 

9.46 

9.99 

10.02 

10.26 

10.36 

A E [GeV] 

0.64 

0.20 

0.05 

1.10 

0.67 

0.54 

0.31 

0.20 

AM [GeV] 

0.02 

-0.03 

0.03 

0.06 

-0.06 

-0.06 

-0.08 

-0.07 

r 0 [fm] 

0.50 

0.72 

0.90 

0.28 

0.44 

0.56 

0.68 

0.78 


Table 3: Quarkonium Spectroscopy from Non-Relativistic Potential Theory (E! 

We thus see that in particular the J/$ and the lower-lying bottomonium states are very 
tightly bound (2Md,b — l ) and of very small spatial size (ro 2r^ ~ 2 fm). 

Through what kind of interaction dynamics can they then be dissociated? 
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As illustration, we consider the collision dissociation of a J/-0. It is very small ( rj ~ 0.25 
fm) and hence can only be resolved by a sufficiently hard probe. It is moreover tightly 
bound (2 M D — Mj/^ ~ 0.6 GeV), so that only a sufficiently energetic projectile can break 
the binding. Hence in a collision with a normal hadron, the J/ip can only be dissociated 
by an interaction with a hard gluon contituent of the hadron, not with the hadron as a 
whole (see Fig. ED. 



Figure 7: J/-0-hadron interaction 


The gluon momentum distribution g(x) in a hadron is determined in deep inelastic lepton- 
hadron scattering; with kh denoting the gluon momentum, x = k^/p^ specifies the fraction 
of the incident hadron momentum p^ carried by the gluon. For mesons, one folds for large 
momenta 

g(x) ~ (1-x) 3 , (14) 

so that the average momentum of a hadronic gluon is 

{k)h = \(Ph)- (15) 

For thermal hadrons in confined matter, (ph) ~ 3 T, with T < 175 MeV, so that with 

{k) h — -T < 0.1 GeV < 0.6 GeV (16) 

5 

the gluon momentum is far too low to allow a dissociation of the J/-0. 

On the other hand, the average momentum of a deconhned thermal gluon in a quark-gluon 
plasma will be 

(k g ) ^ 3 T (17) 

and for T > 1.15 T c , this provides enough energy to overcome the J/'if) binding. We 
thus expect that a hot deconhned medium can lead to J/rp dissociation, while the gluons 
available in a confined medium are too soft to allow this. 

To make these considerations quantitative, one first has to calculate the cross-section for 
the gluon-dissociation of a J/^, a QCD analogue of the photo-effect. This can be carried 
out using the operator product expansion EiEU, and the result is 


1 (k/AE^ — l) 3 / 2 
~ m\ {k/AE^f 


(18) 
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with AEj/^ = 2 M d — Mj/^. The corresponding cross-section for the hadron dissociation 
is then obtained by convoluting the gluon-dissociation cross-section m with the hadronic 
gluon distribution function g(x), which for J/4 >-meson interactions leads to 

<7h-J/1> — CTgeomll ~ A 0 /A) 5 ' 5 (19) 

with A ~ (s — and Ao — ( Mh + A E^). Here cr ge om — 7r~ 2 mb is the 

geometric J/^ cross-section and Mh denotes the mass of the incident meson. In Fig. |H1 
we compare the two dissociation cross-sections m and m as function of the incident 
projectile momentum. 



k [GeV] 

Figure 8: Gluon and hadron J/ip dissociation cross-sections |21| 

This result confirms the qualitative arguments given at the beginning of this section: 
typical thermal gluon momenta near 1 GeV produce a large dissociation cross-section, 
whereas hadron momenta in a thermal range (up to 2 - 3 GeV) still lead to a vanishingly 
small cross-section. In other words, the J/ij) should survive in confining media, but become 
dissociated in a hot quark-gluon plasma. 

The calculations leading to eq. m are exact in the large quark mass limit mg —> oo. It 
is not clear if the charm quark mass really satisfies this condition, and so the dissociation 
cross section in J/tf }-hadron collisions has been discussed in other approaches, some of 
which lead to much larger values |22j . The basic question is apparently whether or not 
the charmonium wave function has enough overlap with that of the usual hadrons to lead 
to a sizeable effect. It may well be that this question is resolvable only experimentally, by 
testing if slow charmonia in normal nuclear matter suffer significant dissociation. Such 
experiments are definitely possible [23]. 

4. Thermal Quarkonium Dissociation 

We shall here address the first of our “basic” questions: how can we identify through 
measurements the states of matter in a fully equilibrated QCD medium? Assume that 
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we are given a box of such strongly interacting matter, in a stationary state of thermal 
equilibrium at a given temperature. Is it possible to specify the state of the given matter 
and transitions between different states in terms of observable quantities calculated in 
QCD? 

4.1 Interaction Range and Colour Screening 

Consider a colour-singlet bound state of a heavy quark Q and its antiquark Q, put into 
the medium in such a way that we can measure the energy of the system as function of 
the QQ separation r (see Fig. 0 . The quarks are assumed to be heavy so that they are 
static and any energy changes indicate changes in the binding energy. We consider first 
the case of vanishing baryon density; at T = 0, the box is therefore empty. 


Q Q Q Q 



Figure 9: String breaking for a QQ system 


In vacuum, i.e., at T — 0, the free energy of the QQ pair is assumed to have the string 
form PH 

F(r) ~ err (20) 

where a ~ 0.16 GeV 2 is the string tension as determined in the spectroscopy of heavy 
quark resonances (charmonium and bottomonium states). Thus F(r) increases with sep¬ 
aration distance; but when it reaches the value of a pair of dressed light quarks (about 
the mass of a p meson), it becomes energetically favorable to produce a qq pair from the 
vacuum, break the string and form two light-heavy mesons (Qq) and (Qq). These can 
now be separated arbitrarily far without changing the energy of the system (Fig. 0. 

The string breaking energy for charm quarks is found to be 

F 0 = 2 (M d — m c ) ~ 1.2 GeV; (21) 

for bottom quarks, one obtains the same value, 

F 0 = 2 (M b - m b ) ~ 1.2 GeV, (22) 

using in both cases the quark mass values obtained in the solution leading to Table 3. 
Hence the onset of string breaking is evidently a property of the vacuum as a medium. It 
occurs when the two heavy quarks are separated by a distance 

r 0 ~ 1.2 GeV/cx ~ 1.5 fm, (23) 
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independent of the mass of the (heavy) quarks connected by the string. 

ff we heat the system to get T > 0, the medium begins to contain light mesons, and 
the large distance QQ potential F(oo,T) decreases, since we can use these light hadrons 
to achieve an earlier string breaking through a kind of flip-flop recoupling of quark con¬ 
stituents resulting in an effective screening of the interquark force (see Fig. HDD- Near 
the deconfinement point, the hadron density increases rapidly, and hence the recoupling 
dissociation becomes much more effective, causing a considerable decrease of F(oo,T). 




Figure 10: In-medium string breaking through recoupling 

A further increase of T will eventually bring the medium to the deconfinement point 
T c , at which chiral symmetry restoration causes a rather abrupt drop of the light quark 
dressing (equivalently, of the constituent quark mass), increasing strongly the density of 
constituents. As a consequence, F(oo,T) now continues to drop sharply. Above T c , light 
quarks and gluons become deconfined colour charges, and this quark-gluon plasma leads to 
a colour screening, which limits the range of the strong interaction. The colour screening 
radius rp, which determines this range, is inversely proportional to the density of charges, 
so that it decreases with increasing temperature. As a result, the QQ interaction becomes 
more and more short-ranged. 

In summary, starting from T — 0, the QQ probe first tests vacuum string breaking, then 
a screening-like dissociation through recoupling of constituent quarks, and finally genuine 
colour screening. In Fig. we show the behaviour obtained in full two-flavour QCD for 
the colour-singlet QQ free energy as a function of r for different T jlZSj. 



Figure 11: The colour singlet QQ free energy F(r,T) vs. r at different T |23. 
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It is evident in Fig. [TT] that the asymptotic value F(oo,T), i.e., the energy needed to 
separate the QQ pair, decreases with increasing temperature, as does the separation 
distance at which “the string breaks”. For the moment we consider the latter to be defined 
by the point beyond which the free energy remains constant within errors, returning in 
section 4.3 to a more precise definition. The behaviour of both quantities is shown in 
Fig. El Deconfinement is thus reflected very clearly in the temperature behaviour of the 
heavy quark potential: both the string breaking energy and the interaction range drop 
sharply around T c . The latter decreases from hadronic size in the confinement region to 
much smaller values in the deconfined medium, where colour screening is operative. 




Figure 12: String breaking potential and interaction range at different temperatures 


The in-medium behaviour of heavy quark bound states thus does serve quite well as probe 
of the state of matter in QCD thermodynamics. We had so far just considered QQ bound 
states in general. Let us now turn to a specific state such as the J/-0. What happens when 
the range of the binding force becomes smaller than the radius of the state? Since the c 
and the c can now no longer see each other, the J/ijj must dissociate for temperatures 
above this point. Hence the dissociation points of the different quarkonium states provide 
a way to measure the temperature of the medium. The effect is illustrated schematically 
in Fig. El showing how with increasing temperature the different charmonium states 
“melt” sequentially as function of their binding strength; the most loosely bound state 
disappears first, the ground state last. 



Figure 13: Charmonium spectra at different temperatures 

Moreover, since finite temperature lattice QCD also provides the temperature dependence 
of the energy density, the melting of the different charmonia or bottomonia can be specified 
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as well in terms of e. In Fig El we illustrate this, combining the the energy density from 
Fig. [Hand the force radii from Fig. El It is evident that although -0' and \c are expected 
to melt around T c , the corresponding dissociation energy densities will presumably be 
quite different. 


r D (T) \/o e /j4 



Figure 14: Charmonium dissociation vs. temperature and energy density 


To make these considerations quantitative, we thus have to find a way to determine the 
in-medium melting points of the different quarkonium states. This problem has been 
addressed in three different approaches: 


• Model the heavy quark potential as function of the temperature, V(r, T ), and solve 
the resulting Schrodinger equation ©• 

• Determine the internal energy U(r,T) of a QQ pair at separation distance r from 
lattice results for the corresponding free energy F(r,T), using the thermodynamic 
relation 


U(r,T ) 


* l d[F(r,T)/T\ >j 


F(r, T)-T 


( dF(r,T) \ 

\ dT )’ 


(24) 


and solve the Schrodinger equation with V (r, T) = U (r, T ) as the binding potential. 


• Calculate the quarkonium spectrum directly in finite temperature lattice QCD. 


Clearly the last is the only model-independent way, and it will in the long run provide 
the decisive determination. However, the direct lattice study of charmonium spectra has 
become possible only quite recently, and so far only in quenched QCD (no dynamical 
light quarks); corresponding studies for bottomonia are still more difficult. Hence much 
of what is known so far is based on Schrodinger equation studies with different model 
inputs. To illustrate the model-dependence of the dissociation parameters, we first cite 
some early work using different models for the temperature dependence of V(r,T), then 
some recent studies based on lattice results for F(r, T), and finally summarize the present 
state of direct lattice calculations of charmonia in finite temperature media. 


4.2 Potential Model Studies 

The first quantitative studies of finite temperature charmonium dissociation were 
based on screening in the form obtained in one-dimensional QED, the so-called Schwinger 
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model. The confining part of the Cornell potential (JHJh V{r) ~ err, is the solution of the 
Laplace equation in one space dimension. In this case, Debye-screening leads to m 

{ 1 _ e ~l*r ) 

V(r,T) ~ err <- , (25) 

l W ) 

where /i(T) denotes the screening mass (inverse Debye radius) for the medium at tem¬ 
perature T. This form reproduces at least qualitatively the convergence to a finite large 
distance value D(oo,T) = cr/n(T), and since fJ,(T) increases with T, it also gives the 
expected decrease of the potential with increasing temperature. Combining this with the 
usual Debye screening for the 1/r part of eq. © then leads to 

f 1 — C~^ r 1 CX O c ^ Q/ 

V(r, T) ~ or <--- - -e^ r = - {l - e~^ r } - -e~» r (26) 

( pr J r p j r 

for the screened Cornell potential. In [221, the screening mass was assumed to have the 
form /i(T) ~ 4 T, as obtained in first lattice estimates of screening in high temperature 
SU(N) gauge theory. Solving the Schrodinger equation with these inputs, one found that 
both the ip' and the Xc become dissociated essentially at T ~ T c , while the J/'ijj persisted 
up to about 1.2 T c . Note that as function of the energy density e ~ T 4 , this meant that 
the J/rj; really survives up to much higher e. 

This approach has two basic shortcomings: 

• The Schwinger form (EHD corresponds to the screening of err in one space dimension; 
the correct result in three space dimensions is different |2ZJ- 

• The screening mass fr(T) is assumed in its high energy form; lattice studies show 
today that its behaviour near T c is quite different 

While the overall behaviour of this approach provides some first insight into the problem, 
quantitative aspects require a more careful treatment. 

When lattice results for the heavy quark free energy as function of the temperature first 
became available, an alternative description appeared |29j . It assumed that in the ther¬ 
modynamic relation (E)) the entropy term —T(dF/dT) could be neglected, thus equating 
binding potential and free energy, 

V(r, T ) = F{r, T) - T(dF/dT) ~ F(r, T ). (27) 

Using this potential in the Schrodinger equation © specifies the temperature dependence 
of the different charmonium masses. On the other hand, the large distance limit of V(r,T ) 
determines the temperature variation of the open charm meson D, 

2Md(T) ~ 2m c + D(oo, T) (28) 

In fig. H21 we compare the resulting open and hidden charm masses as function of tem¬ 
perature. It is seen that the ^ mass falls below 2 Mp around 0.2 T c , that of the x c at 
about 0.8 T c ; hence these states disappear by strong decay at the quoted temperatures. 
Only the ground state J/-0 survives up to T c and perhaps slightly above; the lattice data 
available at the time did not extend above T c , so further predictions were not possible. 
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Figure 15: Temperature dependence of open and hidden charm masses ESI 


The main shortcoming of this approach is quite evident. The neglect of the entropy term 
in the potential reduces V(r,T) and hence the binding. As a result, the D mass drops 
faster with temperature than that of the charmonium states, and it is this effect which 
leads to the early charmonium dissociation. Moreover, in the lattice studies used here, 
only the colour averaged free energy was calculated, which leads to a further reduction of 
the binding force. 

We conclude from these attempts that for a quantitative potential theory study, the 
free energy has to be formulated in the correct three-dimensional screened Cornell form, 
and it then has to be checked against the space- and temperature-dependence of the 
corresponding colour singlet quantity obtained in lattice QCD. 


4.3 Screening Theory 


The modification of the interaction between two charges immersed in a dilute medium of 
charged constituents is provided by Debye-Hiickel theory, which for the Coulomb potential 
in three space dimensions leads to the well-known Debye screening, 


1 

r 



(29) 


where ro = l//r defines the screening radius [3D|- Screening can be evaluated more 
generally EH for a given free energy F(r) ~ r q in d space dimensions, with an arbitrary 
number q. We shall here apply this to the two terms of the Cornell form, with q — 1 for 
the string term, q — — 1 for the gauge term, in d = 3 space dimensions [ 25] , 


We thus assume that the screening effect can be calculated separately for each term, so 
that the screened free energy becomes 


F(r, T) = F s (r, T) + F c (r, T) 


<rrf s {r,T) - -f c (r,T). 
r 


(30) 


The screening functions f s (r,T) and f c (r,T) must satisfy 


f s (r,T ) = f c (r,T) — 1 for T - 0, 

fs(r,T) = f c (r, T) = 1 forr^O, (31) 


since at T = 0 there is no medium, while in the short-distance limit T 1 r 
medium has no effect. The resulting forms are EH 


Fc(r,T) = 


a 


r L 


3 -/w 


+ nr 


0, the 
(32) 
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for the gauge term, and 


Fs(r,T) 


°\ r (I/ 4 ) 

H [2 3 / 2 r(3/4) 


v fW 

2 3 / 4 r (3/4) 


/i 1/4 [(H 2 ] 


(33) 


for the string term. The first term in eq. © gives the constant large distance limit due 
to colour screening; the second provides a Gaussian cut-off in x = fir , since K x / A (x 2 ) ~ 
exp{—x 2 )), which is in contrast to the exponential cut-off given by the Schwinger form 

(EHD- 




(b) 


Figure 16: Screening fits to the QQ free energy F(r,T ) for T > T c (left) and T < T c 
(right) j2Hl 


At temperatures T > T c , when the medium really consists of unbound colour charges, we 
thus expect the free energy F(r,T) to have the form (1HUI) . with the two screened terms 
given by eqs. (J22D and m- I n Fig. HBK , it is seen that the results for the colour singlet 
free energy calculated in two-flavour QCD HU are indeed described very well by this 
form. The only parameter, the screening mass //, is shown in Fig. El and as expected, it 
first increases rapidly in the transition region and then turns into the perturbative form 
H ~ T. 

The behaviour of QQ binding in a plasma of unbound quarks and gluons is thus well 
described by colour screening. Such a description is in fact found to work well also for 
T < T c , when quark recombination leads to an effective screening-like reduction of the 
interaction range, provided one allows higher order contributions in x = fir in the Bessel 
function A 1 / 4 (x 2 ) governing string screening [2%] . The resulting fit to the two-flavour 
colour singlet free energy below T c is shown in Fig. ITfib. using Ky/Qx 2 + (l/2)x 4 ) in eq. 
(1331) : the corresponding values for the screening mass are included in Fig. H71 

With the free energy F(r, T) of a heavy quark-antiquark pair given in terms of the screen¬ 
ing form obtained from Debye-Hiickel theory, the internal energy U(r,T) can now be 
obtained through the thermodynamic relation (1231) . and this then provides the binding 
potential V(r,T) for the temperature dependent version of the Schrodinger equation ®. 
The resulting solution then specifies the temperature-dependence of charmonium binding 
as based on the correct heavy quark potential [TTi] . Let us briefly comment on the ther¬ 
modynamic basis of this approach. The pressure P of a thermodynamic system is given 
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Figure 17: The screening mass fJ,(T) vs. T 23 EEHJ 

by the free energy, P = —F = —U + TS ; it is determined by the kinetic energy TS at 
temperature T and entropy S(T), reduced by the potential energy U(T) between the con¬ 
stituents. In our case, all quantities give the difference between a thermodynamic system 
containing a QQ pair and the corresponding system without such a pair. The potential 
energy of the QQ pair, due both to the attraction of Q and Q and to the modification 
which the pair causes to the internal energy of the other constituents of the medium, is 
therefore given by U. To determine the dissociation points for the different quarkonium 
states, we thus have to solve the Schrodinger equation (JUJ) with V(r,T ) = U(r,T). 

From eqs. (1321) and (1331) we get for the QQ potential 

V(r,T) = V(oo, T) + V(r,T), (34) 

with 

V(oo,T) = ci- - ct/i + T^-[ci^ + a], (35) 

/i al /I 

and Ci = r(l/4)/2 3 / 2 r(3/4), and where V(r,T ) contains the part of the potential which 
vanishes for r —» oo. The behaviour of V(oo,T) as function of the temperature is shown 
in Fig. HH It measures (twice) the energy of the cloud of light quarks and gluons around 
an isolated heavy quark, of an extension determined by the screening radius, relative to 
the energy contained in such a cloud of the same size without a heavy quark. This energy 
difference arises from the interaction of the heavy quark with light quarks and gluons of 
the medium, and from the modification of the interaction between the light constituents 
themselves, caused by the presence of the heavy charge. - The behaviour of V(r,T) 
is shown in Fig. EH f° r three different values of the temperature. It is seen that with 
increasing T, screening reduces the range of the potential. 

The relevant Schrodinger equation now becomes 

{ — V 2 -y(r,T)Wr) - A^(T)$,(r) (36) 

l m c ) 

where 

A Ei(T) = V (oo, T) - Mi - 2 m c (37) 
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Figure 18: The large distance limit of 
the quarkonium potential V(r,T ) |19| 



Figure 19: Variation of V(r,T) with r 
for different T □33 


is the binding energy of charmonium state i at temperature T. When it vanishes, the 
bound state i no longer exists, so that AE t (T) = 0 determines the dissociation temper¬ 
ature Tj for that state. The temperature enters only through the T-dependence of the 
screening mass /i(T), as obtained from the analysis of the lattice results for T(r, T). In 
Fig. [201 we show the resulting binding energy behaviour for the different charmonium 
states, obtained with m c = 1.25 GeV and yfa = 0.445 GeV; in Fig. [2U we show the 
corresponding bound state radii ng. It is seen that in particular the divergence of the 
radii defines quite well the different dissociation points. 
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Figure 20: T-dependence of binding energy for J/-0 (a) and for \ c (b) flTi] 


The same formalism, with mb = 4.65 GeV replacing m c , leads to the bottomonium dis¬ 
sociation points. The combined quarkonium results are listed in Table 4. They agree 
quite well with those obtained in a very similar study based on corresponding free en¬ 
ergies obtained in quenched lattice QCD |J2|, indicating that above T c gluonic effects 
dominate. Using a parametrically generalized screened Coulomb potential obtained from 
lattice QCD results also leads to very compatible results for the Nf = 2 and quenched 
cases pT5l|. We recall here that the main underlying change, which is responsible for the 
much higher dissociation temperatures for the quarkonium ground states, is the use of 
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Figure 21: T-dependence of bound state radii for J/iJj (a) and for Xc ' (b) flH} 


the full internal energy El, including the entropy term, as potential in the Schrodinger 
equation: this makes the binding much stronger. 

We should note, however, that in all such potential studies it is not so clear what binding 
energies of less than a few MeV or bound state radii of several fcrmi can mean in a medium 
whose temperature is above 200 MeV and which leads to screening radii of less than 0.5 
fm. In such a situation, thermal activation Eg can easily dissciate the bound state. 


state 

JMi S) 

Xc(lP) 

^(2 S) 

T(1S) 

Xb(lP) 

T(2 S) 

Xb(2 P) 

T(3S) 

T d /T c 

2.10 

1.16 

1.12 

> 4.0 

1.76 

1.60 

1.19 

1.17 


Table 4: Quarkonium Dissociation Temperatures na 

4.4 Charmonium Correlators 

The direct spectral analysis of charmonia in finite temperature lattice has come within 
reach only in very recent years [2SJ. It is possible now to evaluate the correlation functions 
Gh(t, T) for hadronic quantum number channels H in terms of the Euclidean time r and 
the temperature T. These correlation functions are directly related to the corresponding 
spectral function a h(M,T), which describe the distribution in mass M at temperature 
T for the channel in question. In Fig. [721 schematic results at different temperatures are 
shown for the J/ip and the Xc channels. It is seen that the spectrum for the ground state 
J/ijj remains essentially unchanged even at 1.5 T c . At 3 T c , however, it has disappeared; 
the remaining spectrum is that of the cc continuum of J/ijj quantum numbers at that 
temperature. In contrast, the Xc is already absent at 1.1 T c , with only the corresponding 
continuum present. 

These results are clearly very promising: they show that in a foreseeable future, the 
dissociation parameters of quarkonia can be determined a b initio in lattice QCD. For 
the moment, however, they remain indicative only, since the underlying calculations were 
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Figure 22: J/'ifj and Xc spectral functions at different temperatures 


generally performed in quenched QCD, i.e., without dynamical quark loops. Since such 
loops are crucial in the break-up of quarkonia into light-heavy mesons, final results require 
calculations in full QCD. Some first calculations in two-flavour QCD have just appeared 
pTTTj and support the late dissociation of the J/tp. The widths of the observed spectral 
signals are at present determined by the precision of the lattice calculations; to study 
the actual physical widths, much higher precision is needed. Finally, one has so far only 
first signals at a few selected points; a temperature scan also requires higher performance 
computational facilities. Since the next generation of computers, in the multi-Teraflops 
range, is presently going into operation, the next years should bring the desired results. 
So far, in view of the mentioned uncertainties in both approaches, the results from direct 
lattice studies and those from the potential model calculations of the previous section 
appear quite compatible. 

5. Quarkonium Production in Nuclear Collisions 

The aim of high energy nuclear collisions is to study colour deconfinement and the resulting 
quark-gluon plasma in the laboratory. Here we want to show how quarkonia can be used 
as a probe in this study. The medium to be probed as well as the quarkonium probes are 
produced in the collision, so that we have to address evolution aspects in both cases. After 
a brief section on the evolution of nuclear collisions, we shall first consider quarkonium 
formation in elementary proton-proton collisions and then turn to in-medium phenomena 
of quarkonia in nuclear collisions. 

5.1 Nuclear Collisions 

Starting from the non-equilibrium configuration of two colliding nuclei, the evolution of 
the collision is assumed to have the form shown in Fig. [221 After the collision, there is a 
short pre-equilibrium stage, in which the primary partons of the colliding nuclei interact 
and thermalize to form the quark-gluon plasma. This then expands, cools and hadronizes. 
The partonic constituents in the initial state of the collision are given by the parton distri¬ 
bution functions of the colliding nuclei. To produce a large-scale thermal system, partons 
from different nucleon-nucleon collisions have to undergo multiple interactions. In the cen- 
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Figure 23: Expected evolution of a nuclear collision 


ter of mass of a high energy collision, the incoming nuclei are strongly Lorentz-contracted; 
the resulting parton distribution in the transverse collision plane is schematically illus¬ 
trated in Fig. 1241 The transverse size of the partons is determined by their intrinsic 
transverse momentum, and the number of partons contained in a nucleon is known from 
deep inelastic scattering experiments. The density of partons increases with both A and 
a/s, and at some critical point, parton percolation occurs E3 and global colour connec¬ 
tion sets in. In the resulting connected medium, partons lose their independent existence 
and well-defined origin, so that there is deconhnement, but not yet thermalization. In 
recent years, such partonic connectivity requirements (parton saturation, colour glass con¬ 
densate) have attracted much attention [TT7113K]: they appear to form a prerequisite for 
subsequent thermalization. 



Figure 24: Parton distributions in the transverse plane of a nucleus-nucleus collision 

Assuming that after an equilibration time To a thermalized plasma of deconfined (but 
interacting) quarks and gluons is formed, the initial energy density of the medium can be 
estimated in terms of the emitted hadrons and the initial interaction volume |!J5} 



(38) 


where ( dNf l /dy) y= o specifies the number of hadrons emitted in the nuclear collision per 
unit rapidity at mid-rapidity, and Wh their average energy. The initial transverse size is 
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determined by the nuclear radius Ra , or the nuclear overlap area in non-central collisions. 
The initial longitudinal extension is governed by the equilibration time, which was origi¬ 
nally taken to be about 1 fm ESI- This appears quite reasonable in media, in which the 
intrinsic QCD scale l~ l ~ 1 fm is the relevant measure. For denser partonic systems, one 
can expect a faster thermalization, and so it seems meaningful to use the parton density 
to scale t 0 to higher collision energies. The number of partons emitted by a colliding 
nucleon is about two at SPS energy EH. and thus the parton density in the transverse 
nucleon area is about 1/fm 2 ; hence the above estimate of r 0 appears applicable. The 
number of partons in a nucleon increases with c.rn.s. collision energy roughly as s 1//6 |40j . 
so that the density at RHIC becomes about 2/fm 2 , at the LHC about 6/fm 2 . In Fig. 1251 
we show the initial energy densities for central A — A collisions with A = 200; the range 
shown extends from the value based on To = 1 fm to that obtained if we scale To accord 
to the mentioned parton densities, always with Wh — 0.5 GeV. We note that in all cases 
the energy densities exceed the deconfinement value e(T c ) ~ 0.5 — 1.0 GeV/fm 3 . We also 
indicate which temperatures are expected to be accessible at the different facilities, using 
the results of Fig. |T] 
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Figure 25: Energy density estimates vs. maximum collision energies, for different acceler¬ 
ators, compared to corresponding temperature 

In the following sections, we shall try to show how one can probe deconfinement in nuclear 
collisions, based on the analysis of charmonium production; the extension to bottomonium 
is quite straightforward. First, the production of the probe must be understood when there 
is no bulk medium, so that we begin with charmonium production in elementary collisions. 
Next, one has to understand the modifications arising when the production occurs in a 
confined medium, for which p—A collisions provide the experimental reference. With the 
probe thus prepared and gauged in confined matter, it can be applied in an environment 
in which there might be deconhnement. 

5.2 Hadroproduction of Charmonia 

The hadroproduction of charmonia occurs in three stages. The first stage is the production 
of a cc pair; because of the large quark mass, this process is well described by perturbative 
QCD (Fig. EH). A parton from the projectile interacts with one from the target; the (non- 
perturbative) parton distributions within the hadrons are determined empirically in other 
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Figure 26: Lowest order diagrams for cc production in hadronic collisions, through gluon 
fusion (a,b) and quark-antiquark annihilation (c). 


reactions, e.g. by deep inelastic lepton-hadron scattering. The produced cc pair is in 
general in a colour octet state. In the second stage, it neutralises its colour, leading to 
the third stage, physical resonances, such as J/ij}, Xc or if '. Colour neutralisation occurs 
by interaction with the surrounding colour field; this and the corresponding resonance 
binding are presumably of non-perturbative nature. 

On a fundamental theoretical level, colour neutralization is not yet fully understood. 
However, the colour evaporation model m provides a simple and experimentally well- 
supported phenomenological approach. In the evaporation process, the cc can either 
combine with light quarks to form open charm mesons (D and D) or bind with each 
other to form a charmonium state. The basic quantity in this description is the total sub¬ 
threshold charm cross section S c5 , obtained by integrating the perturbative cc production 
cross section a over the mass interval from 2 m c to 2At high energy, the dominant 
part of S C c comes from gluon fusion (Fig. !2Hhh so that we have 



with g p (x) and gt(x) denoting the gluon densities, x\ and x 2 the fractional momenta 
of gluons from projectile and target, respectively; a is the gg —► cc cross section. In 
pion-nucleon collisions, there are also significant quark-antiquark contributions (Fig.EEb), 
which become dominant at low energies. The basic statement of the colour evaporation 
model is that the production cross section of any charmonium state i is a fixed fraction 
of the subthreshold charm cross section, 


<?i(s) = fi Scc(s), (40) 

where fi is an energy-independent constant to be determined empirically. It follows 
that the energy dependence of the production cross section for any charmonium state 
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is predicted to be that of the perturbatively calculated sub-threshold charm cross section. 
As a further consequence, the production ratios of different charmonium states 


o^Q) 

<L'0) 


A 

fj 


= const. 


(41) 


must be energy-independent. Both these predictions have been compared in detail to 
charmonium and bottomoninm hadroproduction data over a wide range of energies [32J; 
they are found to be well supported, both in the energy dependence of the cross sections 
and in the constancy of the relative species abundances. Let us consider in more detail 
what this tells us about the hadronization of charm quarks. 

We recall that the relative abundances of light hadrons produced in hadron-hadron and 
e + e _ interactions follow the statistical pattern governed by phase space weights mum]: 
the relative production rates are those predicted by an ideal resonance gas at the confine¬ 
ment /deconfinement transition temperature T c ~ 175 MeV. For two hadron species i and 
j that implies at all (high) collision energies 


Ri/j 




(k_ 

dj 



3/2 

exp -{{mi 


m j)/ T c} 


(42) 


for the ratio of the production rates, with di for the degeneracy (spin, isospin) and rn t for 
the mass of species i. For strange particles, the rates (14211 overpredict the experimental 
data; this can, however, be accommodated by a common strangeness suppression factor 
% ~ 0.5, applied as 7" if the produced hadron contains n strange quarks. 

For the hadroproduction of charm, such a statistical description does not work, as seen 
in three typical instances: 

• The total cc cross section increases with energy by about a factor ten between 
a/s = 20 and 40 GeV, while the light hadron multiplicity only grows by about 20%. 
Lienee the ratios for hadrons with and without charm are not energy-independent. 

• From p — p data one finds for J/ip production a weight factor fj/^ ~ 2.5 x 10~ 2 
FT? . Since the subthreshold cc cross section is about half of the single D production 
cross section [43], this implies R(j/^)/d — 10~ 2 ; the ideal resonance gas gives with 
R(j/ 4 >)/d — 1(T 3 a prediction an order of magnitude smaller. Of the total charm 
production, more goes into the hidden charm sector than statistically allowed. 

• For the production ratio of ip' to J/ 1 p, which have the same charm quark infrastruc¬ 
ture, one finds experimentally over a wide range of collision energies /{j/^) — 0.23. 
This energy-independent ip 1 to J/ip ratio can be accounted for in terms of the char¬ 
monium masses and wave functions; it disagrees strongly with the statistical pre¬ 
diction, which gives with -Ry//(_//,/,) ~ 0.045 a very much smaller value. The same 
holds true for the other measured charmonium states, as seen in Table 5, where we 
list experimental results obtained in 300 GeV/c 7r and proton interactions j3Sj. All 
ratios are given relative to directly produced (IS) J/ip states. 

While the first two points could be accommodated by an energy-dependent ‘charm en¬ 
hancement’ factor, the last rules even this out. In the presently available collision energy 
range, charm production in elementary collisions thus does not seem to be of statistical 


26 



nature. It appears to be determined by parton dynamics at an early stage, rather than 
by the phase space size at the confinement temperature. 


ratio 

V>7 ( J /V0 

Xcl/(J/^) 

W( J /7) 

experimental 

0.23 

1.06 

1.50 

statistical 

0.045 

0.113 

0.148 


Table 5: Charmonium production ratios vs. statistical predictions 

Although the colour evaporation model provides a viable phenomenological description 
of the hadroproduction of quarkonia, leading to correct quantitative predictions up to 
the highest energies under consideration, it cannot predict the fractions fi of the hidden 
charm cross sections, and it can even less describe the space-time evolution of colour 
neutralization. For charmonium production in p-A and A-B collisions, the latter is crucial, 
however, and hence a more detailed description of colour neutralization is needed. 



The colour octet model ra proposes that the colour octet cc combines with a soft collinear 
gluon to form a colour singlet ( cc—g ) state. After a rather short relaxation time r 8 , this 
pre-resonant (cc—g) state turns into the physical cc mode by absorbing the accompanying 
gluon, with similar formation processes for Xc and -0' production. The colour octet model 
encounters difficulties if the collinear gluons are treated perturbatively, illustrating once 
more that colour neutralization seems to require non-perturbative elements |4SjJ. However, 
it does provide a conceptual basis for the evolution of the formation process (see Fig. 
EH). The colour neutralization time r 8 of the pre-resonant state can be estimated [4i*)T]: 
it is essentially determined by the lowest momentum possible for confined gluons, r 8 — 
(2m c /) 1//2 ~ 0.25 fm. The resulting scales in J/%1) formation are illustrated in Fig.EBl The 
formation time for the actual physical ground state J/0 is presumably somewhat larger 
than r 8 ; although r j~ r 8 , the heavy c quarks do not move with the velocity of light. For 
the larger higher excited states, the formation times will then be correspondingly larger. 

There is one further important feature to be noted for J/0 hadroproduction. The J/-0’s 
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Figure 28: Scales of J/if; production 


actually measured in hadron-hadron collisions have three distinct origins: about 60 % 
are directly produced IS charmonium states, while about 30 % come from the decay 
X C (1P) —> J/^+ anything, and the remaining 10 % from r//(2 S) —» J/V’+ anything |50j . 
Such feed-down also occurs in T production eh. In all cases, the decay widths of the 
involved higher excited states are extremely small (less than one MeV), so that their life¬ 
times are very long. The presence of any medium in nuclear collisions would therefore 
affect these excited states themselves, and not their decay products. 

5.3 Charmonium Production in p—A Collisions 

We have seen that in the use of charmonia to study nuclear collisions, the creation of 
the medium and the production of the probe lead to two distinct formation scales, r 0 
and Tg. In p — A collisions, the presence of normal nuclear matter can affect charmonium 
production. Here there is no formation time for the medium, so that such collisions 
provide a tool to probe charmonium production, evolution and absorption in confined 
matter. 

Nuclear effects can arise in all the evolution stages of J/-0 production, and a number of 
different phenomena have been studied in considerable detail. 

• The presence of other nucleons in the nucleus can modify the initial state par- 
ton distribution functions, which enter in the perturbative cc production process 
illustrated in Fig. [23 This can lead to a decrease ( shadowing ) or to an increase 
( antishadowing) of the production rate. 

• Once it is produced, the cc pair can suffer absorption in the pre-resonance as well 
as in the resonance stage, caused by the successive interactions with the target 
nucleons. 

• The extent of this absorption can moreover be modified through coherence effects 
for the interactions on different nucleons, which can lead to partial cancellation 
(Landau-Pomeranchuk-Migdal effect). 

In all cases, the crucial quantity is the momentum of the charmonium state as measured 
in the nuclear target rest frame. 

Since we eventually want to probe the effect which the ‘secondary’ medium produced 
by nucleus-nucleus collisions has on charmonium production, it is of course essential to 
account correctly for any effects of the nuclear medium initially present. Let us therefore 
first summarize the main features observed for charmonium production in p —A collision 
experiments [52| . 





• At fixed collision energy, quarkonium production rates per target nucleon decrease 
with increasing A. 

• The production rates decrease for increasing J/-0 momentum as measured in the 
nuclear target rest frame. 

• The nuclear reduction at p — N mid-rapidity appears to become weaker with in¬ 
creasing collision energy. 

• For fixed collision energy, mass number A and J/ip rapidity, the reduction appears 
to increase with the centrality of the collision. 

• At sufficiently high momentum in the target rest frame, the different charmonium 
states appear to suffer the same amount of reduction, while at lower energy, the 
-0' is affected more than the J/'ijj. 

At present, there does not seem to exist a theoretical scenario able to account quantita¬ 
tively for all these observations. In fact, so far not even a common scaling variable for 
the different effects has been found. Shadowing would suggest scaling in the fractional 
target parton momentum x 2 , while absorption of the pre-resonance state would point to 
the fractional beam momentum Xp. Neither is in good accord with the data. For a further 
discussion, we therefore refer to recent reviews IS21ES! and concentrate here on obtaining 
some operational methods to specify nuclear effects in p — A collisions in a way that can 
be extended to A — B collisions. This, incidentally, illustrates the crucial importance of 
having p—A data in order to arrive at an interpretation ot A —A results; for future LHC 
experiments, where parton saturation may play a crucial role, this is even more so the 
case. 

Consider the production of a charmonium state in a p—A collision M- If the initial cc 
production process occurs at a random point inside the nucleus, the evolving charmonium 
state will on the average traverse a path of length La = 3i?^/4 through nuclear matter. 
We have to know in what stage of its evolution the charmonium will be along this path. 
The distance travelled by the charmonium is, in the rest frame of the nucleus, given by 



(43) 


where Pa is the charmonium momentum in the nuclear rest frame, M the mass and r 
the proper age of the charmonium, i.e., its lifet im e in its own rest frame. For simplicity, 
we consider for the moment only charmonia produced at rest in the center of mass of a 
proton-nucleon collision (xp = 0), where most of the data are taken. Then we have 



(44) 


with m for the nucleon mass and a/s for the proton-nucleon c.m.s. collision energy. In 
Table 6, we show as illustration the resulting values for the charmonium momenta Pa in 
the nuclear rest frame at different collision energies; FNAL refers to the energy of the 
fixed target program. We also give the corresponding colour neutralization path lengths 
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d 8 = d{r = r 8 ). Note that these are the values for xp = 0; Pa and d 8 depend quite 
sensitively on xp. 


accelerator 

a/s [GeV] 

Pa [GeV] 

d 8 [fm] 

SPS 

17 

29 

2 

SPS 

29 

48 

4 

FNAL, HERA B 

40 

66 

5 

RHIC 

200 

330 

26 

LHC 

5500 

9000 

730 


Table 6: Charmonium formation parameters at xp = 0 for different collision energies 

We now follow the charmonium in time along its path through the nucleus. For times 
shorter than the colour neutralization time r 8 , r < r 8 , the nucleus is traversed by a pre¬ 
resonance state, which can be dissociated through interaction with nuclear matter, with 
a dissociation cross section Odi S s- For r > r 8 , the physical resonance passes through the 
remaining nuclear medium. From Table 5 we note that for RffIC and LffC, the colour 
neutralization point lies well outside any nucleus, so that there the nuclear medium only 
sees the pre-resonance state. For FNAL and ffERA B, d 8 is approximately equal to the 
size of the largest nuclei, so that here as well only pre-resonance effects should play a role. 
Since up to times r 8 the different charmonium states are indistinguishable, absorption 
effects for yfs > 40 GeV (at xf = 0) should be the same for all states. At the SPS, both 
pre-resonance and resonance absorption can come into play, so that the effect should be 
stronger for the '0 / than for the J/ip. Here it should be noted, however, that in both cases, 
the expansion of the cc occurs with less than the speed of light, so that at r = r 8 the 
colour singlet is not yet of physical size. 

Let us now quantify these considerations SSL According to the Glauber formalism of 
nuclear scattering theory, a cc pair formed at point z 0 in a target nucleus A has a survival 
probability 


Sf — J d 2 bdz pA{b, z) exp j — (A — 1) J dz' pA^b, z') cr l diss ^ , (45) 

where the integration covers the path, at impact parameter b , remaining from Zq out 
of the nucleus, and where cr ^ describes the overall “absorption” effect on the observed 
charmonium state i along the path. The result is then averaged over impact parameter 
and path lengths. The traversed medium of nucleus A is parametrized through a Woods- 
Saxon density distribution Pa{z), and by comparing S) to data for different targets A, 
the effective dissociation cross section can be obtained for the J/i[} and ip' absorption in 
nuclear matter. The effect of the charmonium passage through the nucleus will arise from 


30 











a superposition of the different stages; but if part of the passage is carried out as physical 
resonance, higher excited states should lead to larger absorption cross-sections than the 
much smaller ground state J/-0. In accord with this, the most recent SPS analysis, based 
on p — A collisions with Be, Al, Cu, W and Pb targets, gives 


crj/^, = 4.18 ± 0.35 mb 


(46) 


for the J/'ifb, and 


<jpi = 7.3 ±1.6 mb 


for the ij:' as nuclear absorption cross sections at i/s ~ 29 GeV. 


(47) 


Data from the SPS were taken for a number of different nuclear targets, but so far no 
information on the centrality-dependence of the results has been given. In contrast, RHIC 
experiments have provided only data for d — Au collisions, but as function of centrality 
and for three different rapidity ranges El- In the central rapidity range \y\ < 0.35, the 
data is obtained from e + e~ measurements, in the larger rapidity intervals 1.2 < y < 2.2 
and —2.2 < y < —1.2 through forward and backward p + p~ detection. The results are 
shown in Fig. ESI where 


p , > {da / dy)<i-Au 

d ~ Au[y) ~ N coll {da/dy) p - p 
measures the nuclear modification in the given y range, relative to the corresponding p — p 
cross section scaled by the number N co u of binary nucleon-nucleon collisions. This number 
is obtained from the measured number of participant nucleons through a Glauber analysis 
and is also used to parametrize the centrality dependence. 




4 8 12 16 


Figure 29: J/-0 production in d — Au collisions at RHIC [HZ] 


So far, within the rather limited statistics available, it is not easy to quantify these RHIC 
results enough to obtain a reasonable estimate of normal nuclear effects. As first attempt, 
we adopt a similar parametrization as used for SPS results and apply the well-known 
simplified form of eq. (S3), 

S ~ exp{-n 0 cr d issP}, (49) 


31 













where L denotes the path of the cc in the nuclear medium. Glauber analysis |5SJ provides 
the relation between impact parameter b and the number of collisions N co u, and simple 
geometry gives L 2 = R\ — b 2 in terms of b and the nuclear radius Ra ■ A fit of Eq. (1491) 
to the data of Fig. PHI gives 1 


Giiss (y = 1-8) = 3.1 ± 0.2 mb 
Odiss (V = 0) = 1.2 ± 0.4 mb 

(Tdiss {y = -1.7) = -0.1 ± 0.2 mb (50) 

for the corresponding J/ijj dissociation cross sections; for y = —1.7, there thus are es¬ 
sentially no nuclear modifications. We note that here, as for the SPS case, these cross 
sections are just a global way to account for whatever nuclear effects can arise. A more 
detailed analysis based on shadowing and absorption is given in |~ 59] . 

5.4 Charmonium Production in Nuclear Collisions 

The basic assumption in the attempt to create deconhned matter through nuclear col¬ 
lisions is that the excited vacuum left after the passage of the colliding nuclei forms a 
thermal medium. This picture is schematically illustrated in Fig. OH A charmonium 
state produced in such a collision will in its early stages always be subject to the possible 
effects of the nuclear medium, just as it is in p— A collisions. Knowing the p— A behaviour 
at the corresponding energy is thus a necessary baseline for probing the additional effects 
of the produced medium. 



before collision nuclear medium produced medium 


Figure 30: Collision stages 


The Glauber formalism used above to calculate the survival probability of an evolving 
charmonium state in a p — A collision now has to be extended to A — B interactions 
The survival probability at impact parameter b now becomes 

Sf B (b) = f d 2 s dz dz' pa(s, z)pB(b — s, z') x 

i roo roc 

-(A - 1) / dz A Pa(s , z A ) Vi - (B - 1) / dz B Pb(P - s, z B ) cp > , (51) 

J z Q Jz Q J 

as extension of Eq. (S3). Here Zq specifies the formation point of the ccg within nucleus 
A, Zq its position in B. With the dissociation cross sections a l diss determined in p — A 

1 In the fit, we neglect the most peripheral point at N co u , which corresponds to b > Ra u and is thus 

due to nuclear surface rather than medium effects. 
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collisions, Eqs. ;l Hi lfTTl) . the ‘normal’ survival probability, i.e., that due to only the nuclear 
medium, is fully specihed. 

Since experiments cannot directly measure the impact parameter b, we have to specify how 
Eq. (JSH) can be applied to data. The Glauber formalism allows us to calculate the number 
N^ B (b) of participant (‘wounded’) nucleons for a given collision as function of centrality. 
In A—A collisions, the number of unaffected ‘spectator’ nucleons N A (b) can be determined 
directly through the use of a zero-degree calorimeter, and hence N AA (b) = 2A — 2N A (b) is 
accessible experimentally. For asymmetric AB collisions, the situation is somewhat more 
complex. 

At RHIC energy, we make again use of the simplified form (USD- The geometry connecting 
the impact parameter b and path length L in p — Au and Au — Au collisions is illustrated 
in Fig. |^1 the relation between b and N part is again given by a Glauber analysis jGOj . The 
resulting survival probability is 

S? A (y, JV„) = — - Raa{v,n w ) 

exp{ —no [a diss (y) + a diss (-y)\L} 

corresponding to the fact that, for y ^ 0, the charmonium state passes one nucleus at 
rapidity y and the other at rapidity —y. 




Figure 31: Impact parameter relation between p —A and A —A collisions 


Eqs. (EH) and E2D specify the suppression pattern for charmonium states as predicted 
from normal nuclear effects. To use quarkonia as probes for the produced medium, we 
now have to study how the behaviour observed in A—A collisions differs from this predicted 
pattern. 

Several possible and quite different effects have been considered as consequences of the 
produced medium on quarkonium production. To be specific, we again consider the 
charmonium case for illustration. 

• Suppression by comover collisions: A charmonium state produced in a primary 
nucleon-nucleon collision can be dissociated through interactions with the con¬ 
stituents of any medium subsequently formed in the collision. Such dissociation 
could occur in a conhned m as well as in a deconhned medium . 2T| . 

• Suppression by colour screening: If the produced medium is a hot QGP, it will 
dissociate by colour screening the charmonium states produced in primary nucleon- 
nucleon collisions. Due to the rareness of thermal charm quarks in the medium, the 
separated c and c combine at hadronization with light quarks to form open charm 
mesons [3]. 
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• Enhancement by recombination: In the hadronization of the QGP, charmonium for¬ 
mation can occur by binding of a c with a c from different nucleon-nucleon collisions, 
as well as from the same. If the total number of available cc pairs considerably ex¬ 
ceeds their thermal abundance, statistical recombination enhances hidden relative 
to open charm production, as compared to hadron-hadron collisions 

In addition, the partonic initial state of the colliding nuclei, which leads to the formation 
of the produced medium and to that of charmonium states, will change its nature with 
increasing A and/or y 7 ?; eventually, parton percolation (colour glass formation) can lead 
to a very different medium, with possible effects on production and binding of charnronia 

& 

Is it possible for experiment to distinguish between these different scenarios? Before 
turning to the experimental situation, we want to discuss in some more detail the salient 
features of each approach. 

Suppression by Comover Collisions 

If the charmonium state moves in a random scattering pattern through the produced 
medium, its survival rate is approximately given by 

Si = exp{— (JiiirQ ln[n/7if]} (53) 

with ai denoting the dissociation cross section, n the initial density of the medium after 
a formation time ro, and rif the ‘freeze-out’ density, at which the interactions stop. 

Since the cross section for J/ip break-up through gluon collisions is large [21,] and the gluon 
density high, there will be significant charmonium suppression in a deconfined medium, 
even if this is not thermalized. In an equilibrium QGP, this dissociation is presumably 
accounted for by colour screening, provided the effect of the medium on the width of the 
surviving states is also calculated. 

Charmonium dissociation by interaction with hadronic comovers has received consider¬ 
able attention in the past EH- However, if one restricts the possible densities to values 
appropriate to hadronic matter (n/j<0.5 fm -3 ) and the cross sections to those obtained in 
section 3, the effect of hadron dissociation is negligible. Even a cross section increase to 
the high energy limit still leads to less than 10% effects. More recent analyses [S3] thus 
conclude that a hadronic medium will not result in significant suppression. 

In Fig. E3 we illustrate schematically the overall behaviour expected for J/ij) dissocia¬ 
tion through comover collisions, assuming that beyond a deconfinement threshold, the 
comover density increases with energy density in a monotonic fashion, with little or no 
prior suppression in the hadronic regime. 

Suppression by Colour Screening 

The theoretical basis of this effect has been considered in detail in chapter 5; the colour 
field between the heavy quarks becomes modified due to the presence of a medium of 
unbound colour charges. The results obtained for this effect in statistical QCD are as such 
model-independent, once all calculational constraints are removed. What is speculative 
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Figure 32: J/if) suppression by comover 
collisions 


Figure 33: Sequential J/-0 suppression 
by colour screening 


and model-dependent is its application to nuclear collisions: these do not necessarily 
produce the medium studied in thermal QCD, and the different evolution stages in nuclear 
collisions can introduce factors not present in the study of equilibrium thermodynamics, 
such as the oversaturation of cc pairs just mentioned. 

If the medium produced in high energy nuclear collisions is indeed the quark-gluon plasma 
of statistical QCD, and if charmonium production can be treated as a distinct process 
within such an environment, then the effect of colour screening seems clear. The partition 
of the cc pairs produced in nucleon-nucleon collisions into hidden and open charm is 
non-statistical, favouring the hidden charm sector because of dynamical binding effects. 
Colour screening destroys these and hence moves the partition more and more towards a 
statistical distribution, thus suppressing charmonium production rates relative to those 
observed in elementary interactions. 

A crucial feature of J/-0 suppression by deconfinement is its sequential nature pH, EEJ- In 
the feed-down production of J/^, the produced medium affects the intermediate excited 
states, so that with increasing temperature or energy density, first the J/ip's originating 
from ip 1 decay and then those from y c decay should be dissociated. Only considerably 
higher temperatures would be able to remove the directly produced J/ip’s. Such a step¬ 
wise onset of suppression, with specified threshold temperatures, is perhaps the most 
characteristic feature predicted for charmonium as well as bottom on in m production in 
nuclear collisions. It is illustrated schematically in Fig. [131 where we have defined the 
J/il) production probability to be unity if the production rate suffers only the calculated 
nuclear suppression. 

Enhancement through Recombination 

In charmonium hadroproduction, J/^’s are formed because some of the cc pairs produced 
in a given collision form a corresponding bound state. In a collective medium formed 
through superposition of many nucleon-nucleon (NN) collisions, such as a quark-gluon 
plasma, a c from one NN collision can in principle also bind with a c from another NN 
collision. This ‘exogamous’ charmonium formation can lead to enhanced J/pJ production, 
provided the overall charm density of the medium at hadronization is sufficiently high 
and provided the binding probability between charm quarks from different sources is 
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large enough. 

The production of cc pairs in primary collisions is a hard process and thus grows in A —A 
interactions with the number of nucleon-nucleon collisions; in contrast, the multiplicity 
of light hadrons grows roughly as the number of participant nucleons. Hence the relative 
abundance of charm to non-charm quarks will be higher in A — A than in p — p colli¬ 
sions. Moreover, the cc production cross section increases faster with energy than that 
for light hadron production. The two effects together imply that in the medium produced 
in energetic A — A collisions, the ratio of charm to non-charm quarks is initially much 
higher than in a equilibrated QGP. If this oversaturation of charm is preserved in the 
thermalization process, the combination of random c and c quarks from different primary 
nucleon-nucleon interactions becomes more and more likely with increasing energy. 

Whether or not such an enhancement becomes significant depends on two factors. On 
one hand, the initial charm oversaturation must be preserved, so that the total charm 
abundance in non-thermal. On the other hand, it is necessary that the ‘recombination’ 
of the pairs into charmonia is sufficiently strong. Here it is generally assumed that the 
final hadronization occurs according to the available phase space. Thus the number of 
statistically recombined J/-0’s has the form Nj ~ N^/Nh, growing quadratically in the 
number of produced cc pairs; here Nh denotes the number of light hadrons (or quarks). 
This implies that the hidden to open charm ratio, e.g., Nj/^/Nd ~ N c5 /Nh, increases 
with energy, in contrast to the energy independent form obtained for a fully equilibrated 
QGP, or to the decrease predicted by colour screening. 

Several studies of this effect have led to different noticeable enhancement factors M , in 
the strongest form even predicting an overall enhancement of J/-0 production in A—A 
collisions relative to p —p results scaled by binary collisions. A crucial prediction of the 
approach is the increase of the enhancement with centrality, as shown in Fig. El because 
of the corresponding increase in the number of collisions and hence of the number of cc 
pairs. 



Figure 34: J/'ip enhancement by recombination 

Initial State Suppression 

For two colliding nuclei at large enough A or a/s, the density of partons in the transverse 
plane becomes so large (see Fig. Mi that the partons percolate, producing an inter¬ 
connected network. If the resolution scale of the partons in this network is sufficient to 
resolve charmonia, it could also dissociate them. This scenario is in many ways similar 
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to suppression by colour screening, but it is on the whole more model-dependent. Given 
the average transverse size of the partons in terms of the resolution scale of the parton 
distribution function, it is possible to calculate the percolation point; but its effect on the 
different charmonium states is not a priori evident and can only be addressed in models 

mm- 

In other words, the relation between the onset of charmonium dissociation and the per¬ 
colation point remains so far of phenomenological nature. In contrast to thermal QCD, it 
does not seem possible here to evaluate dissociation points from basic theory. One promis¬ 
ing approach might be the study of quarkonium dissociation in the strong classical colour 
fields created in the percolating medium through the onset of colour glass formation. 

6. Experimental Results from Nucleus-Nucleus Collisions 

Over the past twenty years, great experimental efforts have produced a wealth of data on 
charmonium production in nucleus-nucleus collisions. The pioneering work of the NA38 
and NA50 collaborations at the SPS led to the second generation experiment NA60, which 
is just now announcing its first results. At the same time, RHIC is providing us with a 
first look at how things change when the collision energy is increased by a factor ten. It 
is certainly an exciting time; but although the results may give some indications, more 
information is necessary in order to reach final conclusions. So let us try to summarize 
what we know today and then try to point out what that could mean and what more is 
needed before we can be sure. 

6.1 SPS Results 

Extensive data was taken for S — U collisions (NA38) at a/s = 19.4 GeV and for Pb — Pb 
(NA50) collisions at a/s = 17.3 GeV, for J/ip |2ji and for ip' production HU , as 
well as the mentioned reference data for p — A interactions. The main aspects of the 
NA38/NA50 results are shown in Fig. EH Here the charmonium production rates are 
measured relative to Drell-Yan pair production in the same collision configuration, and the 
resulting ratio is then normalized to the value expected when normal nuclear absorption is 
taken into account, with the cross sections of eqs. (I4(i |147l) . The centrality is experimentally 
determined either by a zero degree calorimeter, which measures the overall spectator 
energy ( Ezdc ) and thus leads directly to the number of participant nucleons, or through 
the transverse energy (E?) or multiplicity of the produced hadrons. The results are shown 
as function of centrality as given by the average length L of nuclear matter traversed, 
which can be determined through a Glauber analysis of either E Z dc or E T data. 

In summary, we note the following behaviour: 

• The rates for J /ip and ip' production in p — A collisions are correctly accounted for 
by normal nuclear absorption. 

• J/ip production in S—U collisions at a/s = 20 GeV is compatible with normal nuclear 
absorption, ip' production in the same S — U collisions, however, shows a further 
reduction. This is the first instance of an ‘anomalous’ suppression, decreasing the 
normal rate by up to a factor five in the most central collisions. 
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Figure 35: J/-0 and i/j' production from 
NA38/NA50 H] 



Figure 36: J/-0 production from 

NA38/NA50/NA51/NA60 |7U] 


• ,//-?/; production in Pb — Pb collisions at yfs ~ 17 GeV as function of centrality 
shows an onset of anomalous suppression. While peripheral results again agree with 
extrapolated p — A results, beyond some centrality there is an additional reduction 
by some 20 - 40 %. The anomalous i^ 1 suppression in Pb—Pb interactions starts 
for more peripheral collisions than that for J/-0, and as function of the available 
medium it is compatible to the S—U results. 

The onset of anomalous J/ij) suppression in Pb — Pb collisions and its absence in S — U 
interactions suggested studying a lighter A—A combination, to obtain further information 
about this phenomenon. The NA60 collaboration has now presented first results on 
J/ij) production in In — In collisions [TU); they are shown in Fig. ES3 together with the 
other J/i)) data taken at the SPS. 

• In In —In collisions there is again an onset of anomalous J/ip suppression, which 
here, however, appears to occur at a centrality at which the corresponding S — U 
data do not show any such effect. The Pb—Pb data is within errors compatible with 
an earlier as well as a later onset position. 

Finally we note that also the transverse momentum behaviour of J/-0 production in p—A 
and A —A collisions was studied at the SPS [H|72j. 

• As shown in Fig. EUl the average (squared) J/ip transverse momentum increases 
with the amount of nuclear matter. A detailed view of the Pb — Pb results in Fig. 
EH suggests a reduced broadening rate for more central collisions. 

6.2 RHIC Results 

In the past months, the PHENIX collaboration at RHIC has presented J/ij) production 
results at y/s = 200 GeV for d — Au , Au — Au and Cu — Cu collisions m- We should 
note that the increase of a factor ten in collision energy corresponds to an increase by a 
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Figure 37: Transverse momentum be¬ 
haviour of J/ij} production :72J 


Figure 38: J// transverse momentum 
behaviour in Pb — Pb collisions 


factor two to four in energy density, depending on the value of r 0 , so that the new data 
smoothly extend the regime studied so far. The data are generally presented in the form 
of a nuclear modification factor Raa, giving the J/tp production rate in Au—Au collisions 
relative to the p—p rate scaled up by the number of binary collisions. 

The results for Au — Au collisions 2 at a/s = 200 GeV are given for two rapidity ranges, 
\y\ < 0.35 and \y\ G [1.2, 2.2], as was the case for the reference d—Au collisions (Fig. I2H1 
The resulting rates Raa are shown in Fig. OH 
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Figure 39: J// production rates at RHIC El) 

The corresponding survival probabilities can now be determined by eq. jHH), making use 
of the nuclear medium modification through the simplified Glauber form (14911 . Since the 
resulting cross sections are so far known with very limited accuracy, we do not include the 
additional error arising from RdAu■ In Fig. OH we show the J// survival probability in 
Au—Au collisions for the two rapidity intervals, as function of the number of participants 
-ZVpart- The suppression is seen to increase with centrality, up to a survival rate of 40 - 
60% remaining for the most central Au — Au interactions. 

2 We shall not consider the results of the Cu — Cu interactions here, since at present we do not have 
access to the corresponding Glauber studies needed for comparison. 
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Figure 40: J/"0 survival probability at Figure 41: J/-0 transverse momentum 

RHIC behaviour at RHIC I7H17H1 


Preliminary RHIC data on the transverse momentum behaviour of J/ip production have 
also been presented ca. As seen in Fig. 0TJ the muon arm results are similar to those 
obtained at the SPS, with increasing broadening from p—p to d—Au and on to Au—Au as 
function of centrality. The e + e _ data have so far very limited statistics and, in particular, 
do not (yet?) show a broadening when going from p—p to d—Au. 

6.3 What does it all mean? 

High energy nuclear collision data thus provide a considerable amount of information on 
the in-medium behaviour of charmonium production. We recall that the raison d’etre for 
the experimental programme was to study in the laboratory the deconfinement transition 
and the new deconhned state of matter predicted by statistical QCD. The crucial question 
in the present context is thus whether the pattern of charmonium suppression observed in 
heavy ion studies is in accord with that calculated in finite temperature lattice QCD. Here 
we should keep in mind that it is not a priori obvious that the medium produced in nuclear 
collisions is in fact the equilibrium QGP studied in statistical QCD. A quantitative relation 
between the charmonium suppression features in nuclear collisions and those calculated 
in lattice QCD would thus provide excellent support for QGP production. So where are 
we now? Does a combination of all the different aspects lead to some first conclusions? 
In this section, we want to address several issues that come up if we try to find a coherent 
interpretation. 

We begin by recalling that if we want to compare data from different nuclear collisions 
and from different collision energies, we have to define some more or less ‘universal’ scale 
variable. We had shown SPS results as function of the effective path length L through 
nuclear matter, since much of the data was presented that way. From an experimental 
point of view, E Z dc data gives the number of participants as the most directly measurable 
quantity, while from a theoretical point of view, density variables, such as participant den¬ 
sity, energy density or parton density, are conceptually preferable. To compare different 
A — B configurations at the same energy, the participant density n part is still meaningful, 
but for different A—B at different collision energies, one has to resort to energy or parton 
densities in order to accommodate the resulting modifications. 

The first question that arises by looking at SPS results concerns the onset of anomalous 
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Figure 43: J/ip production in Pb — Pb , 
S—U collisions ( Et data) and in In—In 
collisions ( Ezdc data) [ZD] 


J/ip suppression in S — U , Pb — Pb and In —In collisions. Do the data as shown in 
Fig. ESI imply different onsets in the different reactions? In Figs. f~H21 and ED we see that 
the discrepancy between S — U and In — In remains, though less pronounced, when we 
determine the centrality through the energy density. We also note, however, that different 
experimental methods of determining the centrality can lead to slight differences. The 
onset of anomalous suppression in the Pb — Pb data agrees very well with that found in 
the In —In data when both are shown in terms of the Ezdc determined centrality (Fig. 
ED, while an apparently later Pb—Pb onset, in accord with the S—U behaviour, is found 
for an i? T -based analysis (Fig. El). Finally we also note that a 10% shift of the S — U 
normalization would remove all differences. Further study here would thus seem very 
appropriate, with the aim of comparing only data in which the centrality determination 
is carried out in the same way (ideally, Ezdc determined centrality) for different collision 
configurations and for both J/ip and ip'. 



1 2 3 4 5 

Figure 44: J/ip suppression as function of energy density 
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In Fig. H3J we show the J/0 survival probability Sj for a combination of SPS and 
RHIC data as function of the energy density, with all data determined through E Z dc 
measurements. For the RHIC data, the energy density values for the different centralities 
are obtained by a Glauber analysis based directly on the PHENIX data, with r 0 = 1 fm 

m- 

The suppression pattern shown in Fig.^Jis fairly consistent with that obtained in chapter 
4 from finite temperature lattice QCD. We had found there a suppression of Xc and 
i[)' just above T c , which implies around e ~ 1 GeV/fm 3 , while the J/0(1S) should survive 
up to e ~ 10 GeV/fm 3 or more. For J/0 production in nuclear collisions, this means 
a suppression of the 40% coming from \c and 0' decay starting around 1 GeV/fm 3 , 
while the 60% directly produced J/0 remain unaffected until much higher e. The onset 
of the anomalous suppression in Fig. |44] occurs indeed at the expected energy density, 
and the J/0 survival probability converges towards some 50 - 60%. Before considering 
the possibility of a next suppression onset at large e, we should recall that the energy 
density values for the RHIC data were obtained with a formation time tq = 1 fm. As 
mentioned, To could well be smaller (parton densities suggest To = 0.5 fm), which would 
move the RHIC points to correspondingly higher e. However, given the present statistics, 
an interpretation of the large e behaviour as onset would in any case seem premature. 



e rGeV/fm 3 1 


Figure 45: Universal J/0 and 0' suppression [77. 


If the interpretation in terms of sequential suppression is correct, the data for 0' suppres¬ 
sion can be used to form 0.4 S0 + 0.6, and this should then coincide with the measured 
J/0 data [ 771 . We show in Fig. E^Jthe SPS J/0 data from Pb—Pb and In—In collisions 
together with the mentioned 0' form. Within the present errors, this universal scaling 
relation seems to hold quite well. Note that here the centrality of the J/0 data is deter¬ 
mined by Ezdc-, that for the 0' by Ex- It is obviously very important to obtain also the 
corresponding 0' data with Ezdc as the centrality variable. 

The sequential suppression pattern expected from statistical QCD means that in most of 
the RHIC data as well as in the central SPS data, what is seen should be the surviving 
directly produced J/0(1S). This also has consequences on the transverse momentum 
behaviour, which thus provide a way to check the conclusion E3- The main effect leading 
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to a broadening of transverse momentum spectra in nuclear interactions is the collision 
broadening of the incident gluons which fuse to make cc pairs. A standard random walk 
analysis gives for the squared transverse momentum of the produced J/ip [7811:79 


{Pr)pA = (Pt)pp + 


(54) 


in p—A and to 

(Pt)aa = ( Pt) pp ± N AA S 0 (55) 

in A—A collisions. Here N A denotes the average number of collisions of a projectile parton 
in the target nucleus A, and N AA the sum of the average number of collisions of a projectile 
parton in the target and vice versa, at the given centrality. The parameter <5 0 specifies 
the average “kick” which the incident parton receives in each subsequent collision. The 
crucial parameters are thus the elementary (Pr) pp from p—p interactions and the value 
of S 0 , determined by corresponding p — A data; both depend on the collision energy. The 
A-dependence of N A as well as the behaviour of N AA as function of centrality can be 
obtained through a Glauber analysis, thus specifying the “normal” centrality dependence 
of ( Pt)aa■ 

In p — A interactions, the average number of collisions of the projectile nucleon in the 
nuclear target can be estimated as dEnn 


N c ~ (3/4) 2 R A n 0 cr, 


(56) 


where a denotes the nucleon-nucleon cross section and the factor (3/4) accounts for the 
(spherical) nuclear geometry. If the parton fusion leading to cc production occurs on the 
average in the center of the target nucleus, N A ~ ( N c — l)/2. However, the presence 
of normal nuclear absorption of the charmonium shifts the fusion point further “down¬ 
stream” , increasing the parton path and hence the number of collisions in the target up 
to a maximum of N A ~ N c — 1. A systematic analysis would require Glauber calculations 
making use of the absorption cross sections measured at the corresponding energies. The 
variation of a and of the charmonium absorption cross sections UHl US] with energy lead 
to an estimate N A ~ 3 ± 1. 

If RHIC as well as central SPS data basically show the surviving direct J/^, then one 
should observe the p? behaviour as given by eqs. © and This again leads to a 

relation between the different data sets from SPS O I3| and RHIC !ZS]- From the RHIC 
p + p~ data, we obtain < p\ > pp = 2.51 ±0.21 GeV 2 and < p\ >dAu= 3.96 ±0.28 GeV 2 ; for 
the latter value, we have taken the average of the positive and negative rapidity ranges. 
With N a ~ 3, we then get 


6 r hic _ 0 5 GeV 2 ( 57 ) 

for the parton broadening at RHIC as shown in the large rapidity data. We note here 
that each kinematic data set requires a separate analysis, which for the muon arm data 
appears possible. The central e + e _ data so far give (with large errors) a ( p which for 
p—p is larger than that for p—A, so that here more statistics are needed. - At SPS energy, 
we have [8U] ( Pt) pp = 1-25 ± 0.05 GeV 2 and ( Pt) p u — 1-49 ± 0.05 GeV 2 . Using again 
N A ~ 3, we get 

5q PS ~ 0.08 GeV 2 (58) 
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for the corresponding parton broadening. With all parameters thus determined, we ob¬ 
tain as a characteristic of the transverse momentum broadening the average number of 
collisions of a projectile and a target nucleon in an AA collisions at the given centrality, 

(K AA ) = {(Pt)aa-(p 2 t ) pp }/^ 0- (59) 

If our interpretation of sequential suppression is correct and the observed (p^) due to 
the unaffected direct J/^’s, then the SPS and the RHIC muon data should follow a 
common curve. We see in Fig. I^fil that this is in fact the case. Once the corresponding 
broadening pattern for the RHIC electron data is determined, it should also follow this 
pattern. - The normal centrality dependence of (N^ A ) has also been calculated directly in 
a Glauber analysis based on the measured normal nuclear suppression CH|; the resulting 
curve is included in Fig. 021 ~ Contributrions from J/ip formation through exogamous 
cc recombination should weaken the increase with centrality for RHIC data compared to 
that from the SPS [STIj . 
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Figure 46: Transverse momentum broadening at SPS and RHIC ra 

We thus find that the overall view of the charmonium data today is more than before in 
accord with the results from statistical QCD. The crucial elements for this are 

• the common onset of anomalous suppression for J/ip and at the dissociation value 
e ~ 1 GeV/fm 3 expected for higher excited charmonium states; and 

• the apparent suppression saturation at about 50 - 70%, following the initial onset, 
suggesting a survival of directly produced J/^’s up to much higher energy densities. 

The first observation revises prior experimental indications for distinct onsets for i// and 
J/ij) suppression, at e ~ 1 GeV/fm 3 and 2.5 GeV/fm 3 , respectively. The second obser¬ 
vation becomes meaningful because of the revision of the theoretical claim of a direct 
J/if) suppression just slightly above T c . 

It seems clear what is further needed to make these indications conclusive. 

• More precise 'ijj' data at SPS energy (and with Ezdc determined centrality) could 
substantiate the scaling behaviour of Fig. 

• More precise J/^ data in the “saturation” regime e > 3 GeV/fm 3 could substantiate 
the direct (IS) survival rate (~ 60%). 
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• An eventual further suppression onset at much higher e (LHC?) could confirm the 
sequential suppression pattern predicted by statistical QCD. 

• Finally, corresponding studies of the bottomonium spectrum with its different sup¬ 
pression onsets could provide decisive substantiation. 

Given the present state of the data, any conclusions must necessarily remain tentative or 
even speculative. How could the picture just presented be experimentally falsified, and 
what would be a possible alternative account? The most frequently considered alternative 
is based on the possibility of J/^ formation by “recombination”, i.e., by the binding of 
c and c quarks originating from different nucleon-nucleon collisions. This approach m 
assumes two distinct mechanisms affecting J/^j production in nuclear collisions JB2]: 

• A charmonium state produced in a given nucleon-nucleon interaction can be disso¬ 
ciated either by comover collisions or by screening effects in the hot medium. This 
leads to the decreasing survival probability which was labelled “thermal dissocia¬ 
tion” in Fig. 031 

• Any two c and c in the resulting medium, whatever their origin, can then at 
hadronization combine to form a charmonium state. An estimate of this recombi¬ 
nation possibility leads in Fig. [2] to the curve labelled “statistical recombination”. 

• Since the charm production rate as a hard process grows faster than the rate for 
light quarks, the recombination probability will increase with centrality as well as 
with collision energy. This will eventually lead in A—A collisions to J/ip production 
rates, relative to overall charm production, which are larger than those in p — p 
collisions. 

At present RHIC energy and centralities, the decreasing rate for “direct” production (i.e., 
in a given nucleon-nucleon collision), due to thermal dissociation, is taken to be just 
compensated by the increasing rate of charmonia formed by combination of charm quarks 
from different nucleon-nucleon collisions, thus leading to an approximately constant “sa¬ 
turation” regime, as seen in Figs. 1401 or l44l 

We should here note two specific aspects inherent in this approach. The assumed form 
for the decreasing thermal dissociation pattern is, as far as we know today (see section 4.) 
not that obtained for a QGP in statistical QCD, with its survival of J/'ip( IS) up to much 
higher temperatures; instead, it is a phenomenological extrapolation of SPS data. As 
mentioned, it is not at all obvious that nuclear collisions do lead to the equilibrated QGP 
studied in lattice QCD, and it is also not obvious that even in such a medium, collision 
effects on a short-time scale could not modify the purely thermal behaviour found there 
(see the remarks at the end of sections 4.3 and 4.4). Nevertheless, the form used at present 
in the recombination approach [H2| has no direct quantitative QCD basis. - Secondly, we 
had seen that the survival probability appears to converge to roughly the fraction (60%) 
of the remaining directly produced J/ip's. This certainly needs to be checked in more 
detail, but, if it were true, it would need to find some explanation in a recombination 
approach. 

The ultimate tests to falsify one or the other approach seem in any case quite clear. The 
most characteristic feature of the recombination approach is an eventual increase of the 
survival probability with centrality. If observed, this would rule out a direct relation be¬ 
tween the experimental pattern and that obtained in statistical QCD. On the other hand, 
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if higher centralities lead to further suppression, and if moreover a similar pattern is found 
for bottomonium production, this would seem to rule out recombination of exogamous 
pairs as significant charmonium formation process. It thus appears that the LHC, with 
its higher energy density access, will still play a decisive role in obtaining the final answer 
to the use of quarkonium spectra as probes of the strongly interacting matter produced 
in nuclear collisions. 

7. Outlook 

The theoretical analysis of the in-medium behaviour of quarkonia has greatly advanced 
in the past decade. Potential model studies based on lattice results for the colour-singlet 
free energy as well as direct lattice calculations appear to be converging, and within a 
few more years, the dissociation temperatures for the different quarkonium states will be 
calculated precisely. Through corresponding calculations of the QCD equation of state, 
these temperatures provide the energy density values at which the dissociation occurs. In 
statistical QCD, quarkonia thus allow a spectral analysis of the quark-gluon plasma. 

The application of this analysis in high energy nuclear collisions is so far less conclusive. 
There exists a wealth of data from different interactions at the SPS, and first RHIC results 
are now also available. These results indicate the production of a hot, dense, strongly 
interacting medium. It is not yet clear, however, to what extent this medium is indeed 
the quark-gluon plasma studied in statistical QCD. Is it possible to identify dissociation 
onsets for different charmonium states? Does the initial (non-thermal) overabundance of 
charm survive a subsequent thermalization? Is the apparent J/^ suppression saturation 
from SPS up to fairly central RHIC collisions the survival of the directly produced J/^’s as 
seen in lattice QCD studies, or is it the production of additional J/^’s through exogamous 
cc pairings? 

It is undoubtedly very challenging to construct models incorporating as many of the 
observed features as possible. Nevertheless it seems worthwhile to note that measurements 
of the relative dissociation points of the different quarkonium states might be our only 
chance to compare quantitative ab initio QCD predictions directly to data. 
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